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The Progenitor Stars of Gamma-Ray Bursts 

S. E. Woosley 1 , A. Heger 2 

ABSTRACT 

Those massive stars that, during their deaths, give rise to gamma-ray bursts 
(GRBs) must be endowed with an unusually large amount of angular momentum 
in their inner regions, one to two orders of magnitude greater than the ones that 
make common pulsars. Yet the inclusion of mass loss and angular momentum 
transport by magnetic torques during the precollapse evolution is known to sap 
the core of the necessary rotation. Here we explore the evolution of very rapidly 
rotating, massive stars, including stripped down helium cores that might result 
from mergers or mass transfer in a binary, and single stars that rotate unusually 
rapidly on the main sequence. For the highest possible rotation rates (about 400 
km s _1 ), a novel sort of evolution is encountered in which single stars mix com- 
pletely on the main sequence, never becoming red giants. Such stars, essentially 
massive "blue stragglers", produce helium-oxygen cores that rotate unusually 
rapidly. Such stars might comprise roughly 1% of all stars above 10 solar masses 
and can, under certain circumstances retain enough angular momentum to make 
GRBs. Because this possibility is very sensitive to mass loss, GRBs will be much 
more probable in regions of low metallicity. 

Subject headings: supernovae, gamma-ray bursts, rotation 



1. Introduction 

During the last seven years, compelling evidence has accumulated that gamma-ray 
bursts (GRBs) of the "long-soft" variety are a consequence of the deaths of massive stars. 
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GRBs occur in galaxies and regions of galaxies where vigorous star formation is going on 
(Vreeswijk et al 2001; Fruchter et al 1999; Fruchter et al, 2005, in preparation, as cited in 
Levan et al 2005). Some, perhaps even a large fraction, are accompanied by supernovae of 
an unusual kind (Galama et al. 1998; Hjorth et al. 2002; Stanek et al. 2003; Zeh et al. 2004; 
Levan et al. 2005). These supernovae have lost their hydrogen envelopes and are observed to 
be Type Ibc. This is consistent with the theoretical expectation that even a relativistic jet 
cannot escape a blue or red supergiant with enough energy to make a GRB that lasts only 
tens of seconds (Woosley 1993; MacFadyen & Woosley 1999; Zhang et al. 2004). It is clear 
therefore that at least some GRBs are made when massive Wolf-Rayet (WR) stars die. It is 
equally clear that not all WR star deaths make GRBs. 

The discriminating characteristics of those stars that do make GRBs are very likely their 
mass and rotation rate. All currently favored models for GRBs require so much rotation that 
it plays a dominant role in the explosion mechanism. This is true of the collapsar model 
(MacFadyen & Woosley 1999) where enough angular momentum must be present to form 
disk around a black hole of several solar masses. It is also true of the "millisecond magnetar" 
model (Usov 1994; Wheeler et al. 2000), where the rotational energy must be sufficient to 
produce magnetic fields of order 10 16 gauss, and of the "supranova" model (Vietri & Stella 
1998) where a neutron star supported, in part, by rotation must exist for an extended 
period. The energies inferred for the "hypernovae" associated with GRBs (Iwamoto et al. 
1998; Woosley et al. 1999; Deng et al. 2005) are ~ 10 52 erg. If the energy is to be derived from 
rotation, the corresponding pulsar, in those models that use pulsars, must have an initial 
rotation rate of £1.5 ms. This is more an order of magnitude faster than the fastest observed 
pulsars and close to what a neutron star can tolerate without deformation or excluding a 
disk. Collapsar models need even more rotation and are favored by high stellar mass, which 
more readily gives the requisite black hole. So, some GRBs are the violent deaths of massive 
WR-stars whose cores are very rapidly rotating. 

This rapid rotation of the inner core is at variance with what is needed to make common 
pulsars with rotation periods £.15 ms (Heger et al. 2005). There are two possibilities. Iron 
cores generally collapse with high rotation rates, but the rotation is damped during or shortly 
after the explosion either by gravitational radiation or by the transport of angular momentum 
to the ejecta. Or the core is already rotating slowly enough that its angular momentum does 
not exceed that of a typical pulsar (~ 5 x 10 47 erg s). The first possibility does not seem 
likely at the present time. Gravitational radiation by the r-mode instability does not slow 
a 1 ms pulsar to a period of 10 ms in a few centuries (Arras et al. 2003) and the dissipation 
of the rotational energy of a 1 ms pulsar by means other than gravity waves or neutrinos 
would give much more energy than is observed in a typical supernova. These are not iron 
clad arguments, but do motivate the study of the second possibility. 
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That too leads to a conundrum. If typical massive star death gives slow pulsars, what 
special circumstances give a GRB? In Heger et al. (2005), we showed that current estimates 
of magnetic torques in the interiors of evolving massive stars (Spruit 2002) led naturally to 
rotation rates of pulsars in the range 10 - 15 ms, just what was needed. But this was for the 
most common variety of supernovae, Type Hp, that result from the deaths of red supergiants 
(RSGs). More massive stars and especially stars that lost their hydrogen envelopes early on, 
had cores that rotated more rapidly. Could it be that some fraction of those massive stars 
that evolve through a giant phase and lose their envelopes either to winds or companion 
stars end up making GRBs. The answer is "maybe", but it's not easy. There are the 
twin problems of magnetic torques and WR-mass loss. Without compelling reasons to the 
contrary, one must employ the same prescription for magnetic torques in the evolution of 
GRB progenitors as for pulsar progenitors, and the core spins down a lot during the RSG 
phase. Second, even if the envelope is removed, the vigorous mass loss of typical WR stars 
still carries away a lot of angular momentum. 

It has been known for some time that if the magnetic torques are negligible, which is 
to say much weaker than estimated by Spruit (2002), it is easy to give GRB progenitors the 
necessary angular momentum (Heger et al. 2000; Hirschi et al. 2004, 2005), but then one 
must invoke another mechanism to slow down typical pulsars. It has also been known that 
if the magnetic torques are included and standard mass loss rates are employed, that it is 
very difficult to make GRBs (Heger & Woosley 2003; Woosley & Heger 2004) from any star 
that either passes through a giant phase or loses appreciable mass as a WR star. 

Here we consider possible resolutions to this dilemma. All of them require rapidly 
rotating stars to begin with and a decrease of up to a factor of 10 in the standard WR-mass 
loss rates currently in use by the community. This may not be so great an adjustment as it 
sounds. Not only does one expect some scatter in the mass loss rates of stars having the same 
mass, composition and angular momentum on the main sequence, but WR stars of lower 
metallicity are known to have lower mass loss rates (Vink & de Koter 2005). Moreover, the 
rate of mass loss may not necessarily equal the (angle-averaged) rate of angular momentum 
loss (Maeder 2002), since the mass loss could occur predominantly at the poles. 

We consider two possibilities, simple helium cores parameterized by their rotational 
speeds (some fraction of break up) at helium ignition and rapidly rotating single stars that 
experience complete mixing on the main sequence. The former set of stars might be il- 
lustrative of the outcome of stellar mergers or other binary activity (Smartt et al. 2002; 
Podsiadlowski et al. 2004; Fryer & Heger 2005; Petrovic et al. 2005). The latter represents, 
so far as we know, a novel suggestion for the evolution of GRB progenitors. 
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2. Models and Physics 

The stars are evolved using the KEPLER implicit hydrodynamics package (Weaver et 
al. 1978) using the same parameters for, and treatment of angular momentum transport and 
mixing as in Heger et al. (2000) and Heger et al. (2005). Additional discussion of the physics 
is given in those references and in Petrovic et al. (2005). For those calculations that include 
magnetic torques, the formalism of Spruit (2002) is employed. 

A major deficiency of the code is that it is one dimensional, so any deformation due 
to rotation is not followed. All rotational quantities such as angular momentum, torques, 
etc. are angle averaged, and the Lagrangian mass shells transport these average quantities. 
The rotation is thus "shellular" . An even greater deficiency is that rotation is treated as a 
passive quantity with no back reaction on the stellar structure, that is the centrifugal term 
is not included in the force equation. This is not a bad approximation so long as the ratio of 
centrifugal force to gravity remains small. In practice, this is true except for the outermost 
layers at the end of helium burning. However, in the most rapidly rotating stars studied 
here, centrifugal force can approach or even mildly exceed unity in a small fraction of the 
mass (§ 3). It is likely that these layers are ejected to form a disk around the star, but our 
present treatment is unable to follow this realistically. 



2.1. Rotation Rate 

Main sequence stars of type O and B are known to be rapid rotators. About 0.3% of 
B stars have rotational speeds on the main sequence in excess of two-thirds times the break 
up speed, and the average is 25% of break up (Abt et al. 2002). The fastest rotating stars 
considered here have about 45% of break up, or rotational speeds on the main sequence in 
excess of 350 km s _1 . Gies & Huang (2005) estimate that this would be several per cent of 
field B stars. This is consistent with our expectation that only a percent or so of stars over 
10 M Q have the special properties required to make a GRB. 

The rotation rate for WR stars is not well determined observationally, but is expected, 
on theoretical grounds, to be rapid, at least for low metallicity (Meynet & Maeder 2005). 
We are interested here in a small fraction of WR stars that may have experienced unusual 
evolution and thus feel justified in assuming large values up to and including those that 
would cause large deformation. 
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2.2. Mass Loss and Angular Momentum Loss 

Our results will be quite sensitive to the rates adopted for mass loss, especially for 
hydrogen- deficient stars. For main sequence stars and red supergiants the mass loss rates 
employed in these studies were taken from Nieuwenhuijzen & de Jager (1990). For WR 
stars, a mass-dependent mass loss rate (Langer 1989) was assumed using the scaling law 
established by Wellstein & Langer (1999), but lowered by a factor 3 (Hamann & Koesterke 
1998) to account for clumping (Nugis et al. 1998), though see Brown et al. (2004). Wind- 
driven mass loss in main sequence stars is believed to be metallicity dependent and a scaling 
law oc \/~Z has also been suggested for hot stars (Kudritzki 2000, 2002; Nugis & Lamers 2000; 
Crowther et al. 2002). Woosley et al. (2002) assumed that the same scaling law holds for 
WR stars (Vanbeveren 2001) and blue and red supergiants as well. "Metallicity" is assumed 
here to be the initial abundance of heavy elements, especially of iron, not the abundances 
of new heavy elements like carbon and oxygen in the atmospheres of WC and WO stars. 
The situation has been recently examined for WR stars of Type WN and WC by Vink & de 
Koter (2005) who find a scaling law oc Z om for metallicities in the range log Z/Z Q = -1 to 
0, with a more gradual decline below -1 for WC stars, but continuing to at least -2 for WN 
stars. Since our stars spend significant fraction of their lifetimes as WN stars and since we 
consider stars only down to metallicity -2, this revised scaling suggests that the winds we 
actually employed in the present study could have been smaller by a factor of ~ Z 1 ^ 3 . This 
would make the production of GRB progenitors more likely, as we shall see. 

It will turn out that we are also sensitive to the mass loss rates for stars that, though still 
on the main sequence because their central hydrogen abundances have not gone to zero, have 
low surface hydrogen abundances due to deep rotational mixing. We have used WR-mass 
loss rates for all stars with surface temperatures over 10,000 K and hydrogen mass fraction 
less than 40%. Such loss rates may well be too large, but the difference in final mass between 
stars where WR mass loss is not implemented until helium ignition is not large 

The most important effect of mass loss in the present context is to carry away angular 
momentum. It is generally assumed that the momentum lost is just the angle-averaged value 
at the surface times the mass loss rate. However, it may well be that precisely those very 
rapidly rotating WR stars we want to consider have quite anisotropic mass loss. The higher 
temperature and luminosity at the poles makes the loss greater there. If this is a significant 
effect, Maeder (2002) has shown that the angular momentum loss might be considerably 
reduced, perhaps even leading to "breakup" in OB stars. 

Here we consider this possibility as just another uncertainty in the mass loss rate which 
might again allow us to use values somewhat smaller than the standard ones. It should 
also be noted that the mass loss rates for GRB progenitors inferred from their afterglows 
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are generally much smaller than the standard values. Soderberg et al. (2005), for example, 
determine a pre-explosive mass loss rate for Type Ic SN 1998bw of 6 x 10~ 6 M yr" 1 , 
even though the accepted mass for the supernova is thought to have been around 10 M 
(Iwamoto et al. 1998; Woosley et al. 1999). The mass loss inferred for Type Ic SN 2002ap 
was 5 x 10~ 7 M yr~ x (Berger et al. 2002). We thus consider models here that, in addition 
to the assumed scaling with metallicity, have mass loss rates reduced by up to a factor of 10 
from their standard values scaled by \/~Z. 

2.3. Composition, Parameters, and the Naming of Models 

For the bare helium cores, a composition of 98.5% He, 1.5% 14 N, and a solar complement 
(Lodders 2003) of elements heavier than neon was adapted. Since we are not interested in 
nucleosynthesis in this study and since the energy generation rate and opacity are essen- 
tially independent of the composition of trace elements throughout most of the star, the 
composition enters in chiefly as a modification of the assumed mass loss rate. The mass 
loss rate in these solar metallicity models was scaled by factors of 1., 0.3, and 0.1 to reflect 
how things might vary with metallicity and the uncertainties discussed in § 2.2. A reduction 
factor of 10 for example, might correspond either to a very metal deficient star of 1% Z 
(assuming mass loss scales as Z - 5 ) or to a star of 10% Z Q that had, for whatever reason, an 
angular momentum loss three times smaller than given by the fitting formula. For helium 
stars, we considered only a single mass, 16 M , which is rather typical, after some mass 
loss, to what has been discussed as a progenitor for the supernovae seen with GRBs. Very 
similar results would characterize helium cores that were within a factor of two of 16. The 
models are defined in Table 1 by the ratio of centrifugal force to gravity at the surface at the 
equator when the star has produced a central carbon mass fraction of 1%. Here that ratio 
is taken to be a large value, 47%, but the ratio decreases to 30% just 0.25 M into the star. 
The surface rotational speed was 800 km s _1 . An additional set of models with the initial 
angular momentum reduced by half, i.e., a surface break up fraction of 26%, gave essentially 
the same presupernova core properties (Table 1). 

For single stars, several choices of mass (12, 16, 35 M ) and metallicity (100%, 10%, 
1% Z ) were explored (Table 2). The name of each model is given by its zero age main 
sequence mass, a letter indicating metallicity ("T" = 1% solar; "O" = 10% solar; and "S" = 
solar), and another letter to distinguish choices for mass loss rate during the WR stage and 
whether magnetic torques were applied. 
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3. Results and Discussion 

Our principal results are given in Tables 1 and 2. For the bare helium cores that did 
not lose a lot of mass (Table 1), the equatorial rotation rate remained a substantial fraction 
of critical throughout the evolution. The maximum specific angular momentum and slowest 
rotational period was always at the surface. This makes it difficult to envision how any binary 
interaction could have imparted a faster rotation to the inner core than what is calculated 
here. 



3.1. Helium Cores 

Helium cores that lost little mass (e.g., Models A, F, G, H, I, N, O, P) remained nearly 
rigidly corotating throughout helium burning and retained a nearly constant angular mo- 
mentum. It is interesting to compare Models HE16I and HE16P which, aside from magnetic 
torques, were otherwise identical models. We focus on the inner 3 M because that is where 
the possibility of a millisecond pulsar or an extreme Kerr black hole will be determined. Be- 
cause there was essentially no mass loss and the coupling between the inner and outer core 
was not strong in Model HE16I, angular momentum was conserved in the inner core through- 
out the evolution. At helium ignition, the angular momentum inside 3 M Q was 2.86 x 10 50 
erg sec; at the presupernova stage it was 2.78 x 10 50 erg sec. The magnetic torques in HE16P 
on the other hand transferred significant angular momentum out of the inner core. Both 
models had the same angular momentum, to within a few percent, at helium depletion (Y c 
= 0.01) and both were rigidly rotating. Contraction to a central temperature of 5 x 10 8 K 
also gave models whose specific angular momenta differed by less than 10% in the inner core. 
However, from carbon ignition (8 x 10 8 K) onwards they diverged. By carbon depletion, the 
model without magnetic torques was rotating three times faster in its inner core and by the 
time the iron core collapsed this factor had become 8. Interestingly, the angular momentum 
of the inner 3 M of the model with magnetic torques decreased by a factor of 2 after carbon 
depletion, that is during the last 0.9 years of the stars life, 

The evolution of the mass shedding helium cores differed in an expectable way. The 
more mass lost, the slower was the rotation of the outer core and, if magnetic torques 
were appreciable, the slower the inner core as well. Consider Model HE16L which included 
magnetic torques, had a mass loss rate 30% that regarded as standard for solar metallicity 
stars, and ended its life with a mass of 9.58 M Q . Almost all the mass loss occurred during 
helium burning. The angular momentum in the inner 3 M of its core was the same as the 
other models at helium ignition (3.32 x 10 50 erg sec) but by helium depletion, it had declined 
by over an order of magnitude to 2.79 x 10 49 erg sec. The entire star was still rotating nearly 
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rigidly. After carbon burning and the accompanying contraction and spin up of the inner 
core, the angular momentum in the same 3 M was reduced by an additional factor of 3 at 
carbon depletion and an additional 30% by the time of iron core collapse (6.86 x 10 48 erg 
sec). At carbon depletion the spread in angular velocity from center to surface varied by a 
factor of 5. 

In general, one sees the tendency of magnetic torques to enforce rigid rotation. This 
extracts angular momentum from the inner core when it contracts and spins up in the post- 
helium burning stages. It also brakes the core when extensive mass loss slows the rotation 
of the outer layers. 

In those models with appreciable mass loss the ratio of centrifugal force to gravity 
decreases with time (Table 1) and is never greater than the initial value. However, for 
models with little or no mass loss and large magnetic torques, the centrifugal forces at the 
surface increase during the latter stages of helium burning and carbon burning to the point 
where they are comparable to gravity. This signals a breakdown in our treatment of rotation, 
but fortunately occurs only in the outer few hundredths (Models O and P) to few tenths 
(Models G and H) of a solar mass. What probably happens here is that the star forms a 
disk. Depending on the viscosity, some of that disk will reaccrete so that one ends up with 
a star, roughly spherical throughout most of its interior, but with highly deformed surface 
layers rotating at break up in the equator. Whether this might have an adverse effect on the 
rotation of the inner core is uncertain, but may have no more effect than simply losing an 
equivalent amount of mass. 



3.2. Rapidly Rotating Single Stars 

In rapidly rotating single stars one has to deal with the additional complexity of a 
hydrogen envelope. In models with typical rotation rates, say 200 km s _1 , for the assumed 
mixing parameters, the star becomes a supergiant sometime during helium burning. For 
stars with very low metallicity, the supergiant may be blue (BSG), but typically it is red 
(RSG) with a radius of several AU. The formation of this extended, high mass, very slowly 
rotating envelope has a great influence on the core, typically braking its rotation rate below 
that required to make a GRB. In fact, the RSG branch of evolution was extensively explored 
by Heger et al. (2005) and, for the given parameterization, was found to give typical periods 
inferred at birth for pulsars. 

The more rapidly rotating models, which are of greatest interest here, bypass giant 
formation by remaining almost completely mixed on the main sequence. Thus, at the end, 
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they resemble the helium cores of the previous section. An interesting exemplary case is 
Model 16TI. This star rotated at about 400 km s -1 on the main sequence and became a WR 
star shortly after central H depletion. Because of its assumed low metallicity, 1% solar and 
mass loss scaling (M oc Z - 5 ), mass loss on the main sequence was negligible and mass loss 
as a WR star was reduced by 10 compared to solar metallicity stars. In this particular model 
the mass loss was decreased by an additional factor of 3 to explore the consequences. The 
final mass of the 16 solar mass star was 13.95 M . It was a WO star with surface abundances 
40% C, 40%O, and 20% He. Most of the presupernova star was predominantly oxygen and 
heavier elements 

Again focusing on the angular momentum in the inner 3 M (see also Fig. 1), the 
angular momentum when the star had burned about half of its hydrogen (X c = 0.40 was 
1.10 x 10 51 erg sec. This declined by about 20% at central hydrogen depletion and by about 
a factor of 4 at helium ignition. This large decrease was caused by core contraction and 
spin up coupled to the outer layers by magnetic torques. However, half way though helium 
burning the angular momentum in the inner 3 M had increased back to 7.2 x 10 50 erg s 
owing to extensive convection redistributing angular momentum in the core. This decreased 
to 4.80 x 10 50 erg sec when the central temperature was 5 x 10 8 K and to 1.38 x 10 50 erg sec 
at carbon depletion. In marked contrast to the RSGs that make pulsars (Heger et al. 2005), 
the angular momentum in the inner core continued to decline appreciably in the late stages 
of evolution. After central carbon depletion, with only 0.16 years left to live, the angular 
momentum in the inner 3 M Q declined an additional factor of almost 4 to 3.67 x 10 49 erg 
sec. at the presupernova stage. In the presupernova (central density 4.2 x 10 9 g cm -3 ), the 
angular velocity in the inner solar mass was 100 times that near the surface. 

To summarize, for very rapidly rotating stars with magnetic torques and little mass 
loss, angular momentum in the inner core is essentially preserved throughout the main 
sequence and has only declines by about a factor of 2 at helium depletion. The major angular 
momentum loss occurs during carbon and oxygen burning with a large fraction occurring 
during the final months of the star's life. This is because of the very large differential rotation 
developed by the core as it contracts through these advanced burning stages and the effect 
of magnetic torques which try to maintain rigid rotation. 

Like the helium cores, some of these rapidly rotating stars also develop centrifugally 
supported surfaces. In the case of Model 16TI, centrifugal forces exceed gravity at the 
equator during helium shell burning, though only in the outer 0.1 M . Other models to 
develop critical rotation after helium burning were Models 12TC, I and J; 16TB, C, and J; 
120C, J, and N; and 160C, J, M, and N. None of the solar metallicity calculations or 35 
M stars developed critical rotation. 
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3.3. Massive Oe and Be Stars? 

Oe and Be stars are a subclass of massive stars that show emission lines, usually taken 
to be indicative of a disk (Hanuschik 1996). There is evidence for rapid rotation. Indeed Be 
stars are the most rapidly rotating of all non-compact stars (Townsend et al. 2004). These 
stars appear may correspond to a phase of spin up caused by mass transfer in a close binary 
system or by an internal redistribution of angular momentum. Not all Be stars are observed 
in binary systems, so it is possible that some form from single B stars. 

The low metallicity, rapidly rotating stars considered here might evolve through a stage 
having properties similar to Oe and Be stars. However, this would only be for stars that 
had, for some reason, very low mass loss rates. It is interesting though that the stars which 
might develop disks are the same stars most likely to produce GRBs. It may be that the 
precursor to a GRB is an Oe or Be star. However, the converse, that all Oe and Be stars 
make GRBs is unlikely, especially in regions with solar metallicity. 

4. Presupernova Characteristics 

The presupernova characteristics of the cores of our models are given in Tables 1 and 2 
and Fig. 2. Various entries give the baryonic mass of the iron core that collapses, its total 
angular momentum, the rotation rate a pulsar would have if the inner 1.7 M collapsed and 
conserved angular momentum, and the Kerr parameter that a black hole would have if it 
formed from the inner 3 M Q of the model. The baryonic mass of the iron core differs from the 
pulsar mass for various reasons (Timmes et al. 1996), especially because of accretion during 
the explosion and neutrino mass loss. The exact relation is unknown because of uncertainties 
in the explosion mechanism. However, a 1.7 M (baryonic mass) core would give a 1.44 M e 
(gravitational mass) neutron star after neutrino losses. Approximately 20% could be added 
to the rotational period of those models that give neutron stars and not black holes because 
of the angular momentum carried away by the neutrinos (Heger et al. 2005). 

Stars that have an entry greater than 1 for the Kerr parameter at 3 solar masses would 
have to form a disk to carry the extra angular momentum and are thus good candidates for 
collapsars. If a black hole forms in these systems, so will a disk. Other models having a ^0.3 
at 3 M are also good candidates because the angular momentum increases outwards. Fig. 
2 shows the angular momentum distribution in Models 16TI and 160M. Both would form 
accretion disks at 3.5 and 5.5 M Q respectively, even though the Kerr parameter at 3 M Q is 
only 0.44 and 0.25. 
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5. Implications for Gamma-Ray Bursts and Supernovae 

One can broadly characterize the effect of rotation on the explosion mechanism by 
the rotational energy the resulting pulsar would have if one formed and conserved angular 
momentum. For a typical neutron star radius (12 km) and gravitational mass (1.4 M ), 
the moment of inertia is « 0.35MR 2 = 1.4 x 10 45 g cm 2 (Lattimer & Prakash 2001). The 
rotational energy, l/2Iu 2 , is then E rot m 1.1 x 10 51 (5 ms/P) 2 erg. Since the typical kinetic 
energy of a supernova is 10 51 erg, this implies a necessary condition that the pulsar contribute 
the bulk of the energy is that its rotation rate be £5 ms. This lays aside all considerations of 
how this energy might be tapped and with what efficiency To give a "hypernova" with ten 
times this energy requires rotational periods ^2 ms. It turns out that stars that would give 
disks around black holes of several solar masses also require a comparable rotation rate, si 
ms, though the relevant angular momentum is located somewhat farther out in the star. So 
we can make the distinction. Neutron stars with initial periods of 10 ms or longer probably 
won't have a large effect on the explosion; those with periods less than 5 ms might, and a 1 
ms period is needed to make a GRB. 

By this criterion, all neutron stars resulting from stars that pass through a supergiant 
phase (either red or blue) and do not lose their envelopes will rotate too slowly to be GRBs, or 
even to power normal supernovae. Even those stars resulting from helium cores rotating near 
break up - whether formed from binary evolution, or very rapidly rotating solitary stars - will 
be too slow to make GRBs unless their mass loss rates are smaller than generally assumed. 
This could come about either because the metallicity of GRBs is quite low or because the 
rates currently in use overestimate the actual angular momentum loss for unknown reasons 
(§ 2.2. The recent revision downwards of mass loss rates for metal deficient WR stars by 
Vink & de Koter (2005) is helpful in this regard. 

With reasonable variations then it is possible to produce a subset of models that do give 
GRBs. This could be the small number of O and B stars out on the tail of the rotational 
velocity distribution with v rot ~ 400 km s" 1 , which might be a few percent of all such 
stars. It could equally well be a population of helium cores in binaries that have arrived at 
helium ignition with a rotation rate corresponding to about one-third of break up - provided 
the statistics yield a sufficient number of such objects. Either way, because of the likely 
dependence of mass loss on metallicity (Vink & de Koter 2005), GRBs will be favored in 
regions of low metallicity as predicted by MacFadyen & Woosley (1999) and observed by 
Prochaska et al. (2004). 
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Table 1. Presupernova Models For Bare Helium Cores 



Mass/ 


J initial 


Percent 


M 


B-field 


M fina i 


Fe-core 


J Fe — core 


Period 


a BH 


Model 


(10 52 erg s) 


break-up 8. 


WR 




(M ) 


(M ) 


(10 47 erg s) 


(ms) 


(3 M ) 


HE16A 


2.0 


50 


0. 


no 


15.70 


1.91 


1870 


0.06 


(6.0) 


HE16B 


2.0 


7 


1.0 


no 


5.10 


1.45 


275 


0.27 


(1.0) 


HE16C 


2.0 


4 


1.0 


yes 


5.15 


1.44 


5.8 


12.0 


0.02 


HE16D 


2.0 


18 


0.3 


yes 


9.53 


1.61 


33.4 


2.4 


0.16 


HE16E 


2.0 


33 


0.1 


yes 


12.86 


1.38 


45.4 


1.4 


0.23 


HE16F 


2.0 


42 


0.03 


yes 


14.80 


1.75 


114 


0.86 


0.48 


HE16G 


2.0 


47 


0.01 


yes 


15.56 


1.90 


146 


0.82 


0.51 


HE16H 


2.0 


50 


0. 


yes 


15.68 


1.92 


149 


0.82 


0.51 


HE16I 


1.0 


26 


0. 


no 


15.88 


1.89 


1060 


0.11 


(3.6) 


HE16J 


1.0 


5 


1.0 


no 


5.13 


1.42 


245 


0.29 


0.91 


HE16K 


1.0 


2 


1.0 


yes 


5.16 


1.44 


5.7 


12.0 


0.02 


HE16L 


1.0 


10 


0.3 


yes 


9.58 


1.62 


23.6 


3.5 


0.11 


HE16M 


1.0 


17 


0.1 


yes 


13.04 


1.45 


44.4 


1.6 


0.19 


HE16N 


1.0 


22 


0.03 


yes 


14.95 


1.54 


69.0 


1.0 


0.34 


HE160 


1.0 


21 


0.01 


yes 


15.62 


1.64 


87.4 


0.97 


0.42 


HE16P 


1.0 


26 


0. 


yes 


15.88 


1.83 


119 


0.91 


0.45 



a All models had a surface rotation rate of 47% (Models A-H) or 26% (Models I-P) critical 
at helium ignition. The value given here is the percentage of critical rotation at the surface 
when the star has burned half its helium at its center. 
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Table 2. Presupernova Models for Rapidly Rotating Stars 



Mass/ 


T. . a 

° mit 


v rot 


PreSN 


M 


B-field 


Mfi na i 


Fe-core 


JFe— core 


Period 


aBH 


Model (10 


32 erg s) 


(km s _1 ) 




WR 




(M ) 


(M©) 


(10 47 erg s) 


(ms) 


(3 M 


12TA 


0. 


0. 


RSG 




no 


11.96 


1.35 


0. 


_ 


_ 


12TB 


1.5 


290 


RSG 




no 


11.95 


1.42 


362 


0.17 


(2.4) 


i OTP 


z.u 


380 

OoU 


VV It 


n i 

U. 1 


no 


11 

11.00 


1 K7 
1.0 ( 


1 1 on 
i iyu 


n n7 

U.U ( 


f3 7^ 


1 9TD 

1Z 1 YJ 


n 


n 
u. 


rto^or 




no 


11 OK 

i i.yu 


1 3^ 
1 .oo 


n 
u. 






1 9TT? 
iz i ri; 


1.0 


900 

zyu 


ItDVj 




no 


11 QK 

ii. yo 


1 .40 


3fiQ 

ooy 


17 
U.l 1 


(0 A\ 


1 9TT? 


o n 
z.u 


380 

OoU 


VV It 


1 n 

l.U 


no 


o 1 a 
y . io 


1 c;n 
1 .ou 


ouo 


O 1 3 
U. 10 


(0 1 


i OTP 


1.0 


900 

zyu 


ItOVj 




yes 


1 1 OR 

ii. yo 


i /i 1 

1 .41 


P. 9 

o.z 


Q K 

y .o 


O OR 

u.uo 


1 9TU 

izin 


9 n 
z.u 


380 

OoU 


VV It 


1 n 

l.U 


yes 


Q 93 

y .zo 


1/18 
1 .4o 


/I 3 9 
40. Z 


1 8 
1 .0 


17 
U. 1 ( 


1 OTT 
IZ 1 1 


o n 
z.u 


380 
OoU 


AAfR 
VV It 


n 3 
u.o 


yes 


i n 7Q 
iu. ( y 


1 K7 
1.0 i 


fit; n 
00. u 


1 9 

1 .z 


O 31 
U.01 


1 9T T 
IZ 1 J 


9 n 
z.u 


380 

OoU 


VV It 


n i 

U. 1 


yes 


1 1 .041: 


1 89 
1 .oZ 


1 

10U 


O 79 
U. ( z 


U.O ( 


16TA 


0. 


0. 


BSG 




no 


15.95 


1.44 


0. 


- 


- 


16TB 


2.5 


305 


WR 


0.1 


no 


15.29 


1.75 


1640 


0.06 


(5.7) 


1 KTP 
10 1 O 


3 3 
0.0 


300 
oyu 


WR 
VV It 


n i 

U. 1 


no 


1 ^ 93 
10. ZO 


1 R1 
1.01 


1 /I 90 
14ZU 


O OR 

u.uo 


(P. C\\ 

(b.U) 


10 1 D 


n 
u. 


n 
u. 


R.QP 
Dob 




no 


1 c; o^ 
io. yo 


1/1/1 
1 .44 


n 
u. 






10 1 Hi 


9 

z.o 


30^ 
ouo 


AA/R 
VV It 


i n 

l.U 


no 


ii no 

ii. yo 


1 A7 
1 .4 I 


A 

4oy 


17 
U.l 1 


(1 A\ 
(1.4; 


10 1 r 


3 3 
O.O 


300 
oyu 


AAfR 
VVrt 


1 n 

l.U 


no 


1 1 37 
11.0 ( 


1 8 K 

l.oO 


oyo 


O 1 8 

U.lo 


(l.b) 


i «TP 

10 1 Vor 


9 

Z.O 


30^ 
ouo 


DOp 




yes 


1 K P.7 
10.0 ( 


1 7^ 
1.(0 


1 Q 9 

iy.z 


K 1 
O.l 


o oc; 
u.uo 


i «TTT 
io i n 


3 3 
0.0 


3Q0 

oyu 


AA/R 
VV It 


1 n 

l.U 


yes 


1 1 .41:0 


1 81 
l.ol 


P.A 1 
04. 1 


1 7 
1. 1 


O 9/1 
U.Z4 


1 RTT 
1011 


3 3 
0.0 


300 
oyu 


AA/R 
VV It 


n 3 
u.o 


yes 


1Q QK 

io. yo 


1 P.r\ 
1 .ou 


8£ 7 
oO. ( 


O QO 

u.yu 


O A A 
U.44 


1 KT T 
10 1 J 


3 3 
0.0 


300 
oyu 


AATR 
VVrt 


n 1 

U. 1 


yes 


1 CC. 91 

10. zl 


1 88 
l.OO 


1 78 
I/O 


O K7 
U.O t 


O R1 
U.01 


120A 


0. 


0. 


RSG 




no 


11.86 


1.25 


0. 






120B 


1.5 


245 


RSG 




no 


11.83 


1.43 


345 


0.18 


(2.2) 


120C 


2.0 


325 


WR 


0.1 


no 


10.98 


1.57 


1120 


0.07 


(5.2) 


120D 


0. 


0. 


RSG 




no 


11.86 


1.32 


0. 






120E 


1.5 


245 


RSG 




no 


11.83 


1.43 


343 


0.18 


(2.2) 


120F 


2.0 


325 


WR 


1.0 


no 


7.61 


1.42 


354 


0.18 


(1.2) 


120G 


1.5 


245 


RSG 




yes 


11.83 


1.41 


6.1 


9.7 


0.06 


120H 


2.0 


325 


WR 


1.0 


yes 


7.69 


1.39 


10.0 


6.1 


0.04 
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Table 2 — Continued 



Mass/ 


T. . a 

•J mit 


V rot 


PreSN 


M 


B-field 


M fina l 


Fe-core 


J Fe— core 


Period 




Model (10 


3 erg s) 


(km s ) 




WR 




(M ) 


(M ) 


(10 4 ' erg s) 


(ms) 


(3 M 


12DT 

X Zj V ./ 1 


2 D 

Zj . W 


325 


WR 


3 

w . u 




Q 71 


1 41 

X X 


30 ^ 


1 6 

X . \J 


1 7 

U.l 1 






4DD 


WR 


n i 

U. 1 


1 1 f~l 




1 ^9 

1 . UZj 


1 1 6D 


n D7 




1 9DK 

x zj yy iv 


Zj . U 


Ann 


WR 

VV IV 


1 n 

X . U 


LIKJ 


7 SO 


1 ^ 

1 .oo 


SI 4 


n 1 h 




120L 


2 5 

Zj . O 


400 

TI V/ VI 


WR 


1 

X . VI 




7 35 


1 35 


q 5 


5 8 


04 


1 90M 

x zj vy iv j. 


2 5 

zj . 


400 

t:UU 


WR 

VV 11 


3 

w . u 




q 50 

Zj . U W 


1 86 

X . (J VJ 


1 00 


1 

X . w 


30 


12QN 

X Z. \7 ^ N 


2 5 

zj . 


400 

t:UU 


WR 

V V X V- 


1 

W. X 




1 03 


1 ^0 


74 1 

I t: . X 


1 

X . w 


44 


1 6DA 

X U V / i X- 


n 

U . 


n 

U . 


RSG 

X V kJ V ' 






1 5 80 


1 43 


n 

V7 . 






1 6DR 


2 5 

zj . xj 


255 


RSG 

X VkJ VI 




n n 


1 5 57 

XtJ . tj ( 


1 84 


^93 

?JZj O 


20 

W . Zj W 




160C 


3 3 

U . U 


325 


WR 

VV 11 


1 

U . X 




14 96 

XT: . Zj VI 


1 63 


1 360 


06 


f4 41 


1 6DD 

X VI \ ' \.J 


n 

U . 


n 

U . 


RSG 

X VkJ VI 




n n 


1 5 80 


1 43 

X .10 


n 






160E 


2.5 


255 


RSG 




no 


15.57 


1.84 


523 


0.20 


(1.8) 


160F 


3.3 


325 


WR 


1.0 


no 


8.97 


1.35 


318 


0.18 


(1.1) 


1 600 


9 5 

Zj . U 


955 

ZjUU 


RSO 

1VO vj 




yes 


1 d.UU 


i 


q 6 


7 n 


u .UU 


16DH 

1 UWll 


3 3 

U . U 


39 s ) 


WR 

V V X V- 


1 

X . w 


yeb 


q is 

C/ . X (J 


1 4^ 


q 8 


7 o 


03 

u . UU 


1 6DT 

X U V / 1 


3 3 

U . U 


39 5 

■ J Zj ' J 


WR 

VV 11 


3 

w . u 


jeb 


19 91 

X Zj . Zj X 


1 6^) 


^^i 3 


1 ^ 

X . o 


26 

U . ZjVJ 


1 fif) T 

1UVJJ 


4 1 

■4:. 1 


'rfcUU 


VV IV 


n i 




±■4:. ZjU 


1 ^fi 


1 9Qf) 


n nfi 

u . uu 


l°- u J 


160K 

X V .' V / 1 V 


4.1 


400 

t: V/ VJ 


WR 

VV 11 


1 

X . VI 


11U 


8 58 


1.52 


3qq 


0.21 




16DL 


4 1 

t: . X 


400 


WR 

V V 1 V 


1 

X . vi 




8 68 


1 ^2 

X . \J Zj 


1 4 q 

It:. J 


^ 6 


0^1 

u . uu 


1 ROM 

X \j vy ivx 


4 1 

t: . X 


400 

t:UU 


WR 

V V 1 V 


3 

w . u 




1 1 q4 

X X . i./ t: 


1 ^ 


^3 3 


1 4 

X .1 


95 

U . zjtj 


1 RON 

1U WIN 


4 1 

4:. X 


zinn 


VV IV 


n i 

U. 1 


yes 


14 17 


1.IO 


191 

IZj ± 




n 

U.^diU 


12SA 


0. 


0. 


RSG 




no 


10.88 


1.40 


0. 






12SB 


2.0 


280 


WR 


0.1 


no 


q.3q 


1.40 


410 


0.16 


(1.7) 


12SC 


3.0 


405 


WR 


0.1 


no 


8.q2 


1.66 


751 


0.11 


(2.8) 


12SD 


0. 


0. 


RSG 




no 


10.88 


1.40 


0. 






12SE 


2.0 


280 


WR 


1.0 


no 


7.27 


1.42 


355 


0.10 


(1.1) 


12SF 


3.0 


405 


WR 


1.0 


no 


4.06 


1.51 


333 


0.23 


(1.1) 


12SG 


2.0 


280 


BSG 




yes 


7.57 


1.57 


0.4 


8.3 


0.04 


12SH 


3.0 


405 


WR 


1.0 


yes 


5.43 


1.46 


6.6 


11 


0.03 
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Table 2 — Continued 



Mass/ 


J init 


v rot 


TT» OAT 

PreSN 


M 


B-neld 


TV jv 

Mfi na l 


re- core 


T 

" Fe— core 


reriod 


aBH 


Model (10 


erg s) 


(km s ) 




WR 




/TV X \ 

( M o) 


( M 0) 


(10 erg s) 


(ms) 


(3 M 


12SI 


3.0 


405 


WR 


0.3 


ves 


6.95 


1.39 


7.1 


8.9 


0.03 


12SJ 


3.0 


405 


WR 


0.1 


ves 
y 130 


9.03 


1.46 


36.6 


1.7 


0.15 


16SA 


0. 


0. 


RSG 




no 


14.65 


1.38 


0. 






16SB 


2.5 


215 


RSG 




no 


14.03 


1.44 


332 


0.23 


(1.2) 


16SC 


3.3 


270 


WR 


0.1 


no 


11.58 


1.73 


492 


0.19 


(1.4) 


16SD 


0. 




RSG 




no 


14.65 


1.38 


0. 






16SE 


2.5 


215 


RSG 




no 


14.03 


1.44 


332 


0.23 


(1 2) 


16SF 


3.3 


270 


WR 


1.0 


no 


6.25 


1.35 


281 


0.20 


(1 1) 


16SG 


2.5 


215 


RSG 




yes 


11.98 


1.49 


9.5 


7.9 


0.04 


16SH 


3.3 


270 


WR 


1.0 


yes 


7.70 


1.35 


8.2 


6.7 


0.03 


16SI 


3.3 


270 


WR 


0.3 


yes 


9.85 


1.72 


14.8 


6.4 


0.06 


16SI 


4.5 


360 


WR 


0.1 


no 


11.10 


1.76 


550 


0.18 


(1.6) 


16SJ 


4.5 


360 


WR 


1.0 


no 


5.33 


1.46 


293 


0.25 


(1.0) 


16SK 


4.5 


360 


WR 


1.0 


yes 


6.30 


1.33 


5.9 


9.7 


0.02 


16SL 


4.5 


360 


WR 


0.3 


yes 


8.31 


1.52 


10.9 


7.6 


0.04 


16SM 


4.5 


360 


WR 


0.1 


yes 


11.22 


1.90 


50.4 


2.3 


0.16 


350A 


14. 


380 


WR 


1.0 


yes 


12.86 


1.52 


13.2 


6.0 


0.05 


350B 


14. 


380 


WR 


0.3 


yes 


21.24 


2.05 


151 


0.86 


0.38 


350C 


14. 


380 


WR 


0.1 


yes 


28.07 


2.02 


230 


0.59 


0.53 


350A 


14. 


380 


WR 





yes 


34.4 


1.29 


193 


0.28 


(1.2) 



Total angular momentum of the zero age main sequence star. 
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Fig. 1. — Distribution of the specific angular momentum with mass in four 16 solar mass 
models evaluated at the zero age main sequence (top line, solid); central helium depletion 
(dotted line, next down); carbon depletion (dashed line, second down); and at the presuper- 
nova stage (dark solid line, bottommost). The upper two models had reduced metallicity 
(1% solar for the left frame; 10% for the right) and mass loss during the WR stage (30%). 
Both rotated at about 400 km s" 1 on the main sequence and avoided forming a red giant. 
The models in the bottom two panels had solar metallicity and more moderate rotation on 
the main sequence. The one on the right rotated faster and ended up a WR star. The one on 
the left died as a red supergiant. The pulsars derived from the top two models rotated 5 to 
10 times faster. Avoiding the formation of a giant star and suppressing mass loss amplifies 
the final rotational momentum of the core. 
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Model 16TI 



Pre-Supernova 
Angular Momentum 

J I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



2 4 6 8 10 12 14 

m / Mq 




Pre-Supernova 
Angular Momentum 

J I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



2 4 6 8 10 12 14 

m / Mq 

Fig. 2. — Angular momentum distribution of Models 16TI and 160M at the time the star 
collapsed. The dark solid line shows the distribution of specific angular momentum in the 
presupernova star. The blue line indicates the angular momentum required to support matter 
at the stable orbit for a black hole that is not rotating; the green line is for a Kerr black hole 
with rotational parameter a = 1. The red line indicates the last stable orbit for a black hole 
with the mass and angular momentum inside the indicated coordinate in the presupernova 
star. Where the black line is above the red line a disk can form and collapsars are possible. 



